BREATHING MODES OF LONG JOSEPHSON JUNCTIONS WITH 

PHASE-SHIFTS* 
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Abstract. We consider a spatially inhomogeneous sine-Gordon equation with a time-periodic 
drive, modeling a microwave driven long Josephson junction with phase-shifts. Under appropriate 
conditions, Josephson junctions with phase-shifts can have a spatially nonuniform ground state. In 
recent reports [Phys. Rev. Lett. 98, 117006 (2007), arXiv:0903.1046 , it is experimentally shown 
that a microwave drive can be used to measure the eigenfrequency of a junction's ground state. 
Such a microwave spectroscopy is based on the observation that when the frequency of the applied 
microwave is in the vicinity of the natural frequency of the ground state, the junction can switch to 
a resistive state, characterized by a non-zero junction voltage. It was conjectured that the process 
is analogous to the resonant phenomenon in a simple pendulum motion driven by a time periodic 
external force. In the case of long junctions with phase-shifts, it would be a resonance between 
the internal breathing mode of the ground state and the microwave field. Nonetheless, it was also 
reported that the microwave power needed to switch the junction into a resistive state depends on 
the magnitude of the eigenfrequency to be measured. Using multiple scale expansions, we show 
here that an infinitely long Josephson junction with phase-shifts cannot be switched to a resistive 
state by microwave field with frequency close to the system's eigenfrequency, provided that the 
applied microwave amplitude is small enough, which confirms the experimental observations. It is 
because higher harmonics with frequencies in the continuous spectrum are excited, in the form of 
continuous wave radiation. The breathing mode thus experiences radiative damping. In the absence 
of driving, the breathing mode decays at rates of at most 0{t~^^'^) and for junctions with 

a uniform and nonuniform ground state, respectively. The presence of applied microwaves balances 
the nonlinear damping, creating a stable breather mode oscillation. As a particular example, we 
consider the so-called — vr — and — k Josephson junctions, respectively representing the two 
cases. We confirm our analytical results numerically. Using our numerical computations, we also 
show that there is a critical microwave amplitude at which the junction switches to the resistive 
state. Yet, it appears that the switching process is not necessarily caused by the breathing mode. 
We show a case where a junction switches to a resistive state due to the continuous wave background 
becoming modulationally unstable. 

Key words, radiative damping, fractional fiuxon, breathing modes, Josephson junction with 
phase-shifts, long Josephson junction. 
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1. Introduction. A Josephson junction is an electronic circuit consisting of two 
superconductors connected by a thin nonsuperconducting layer, and is the basis of 
a large number of developments both in fundamental research and in application 
to electronic devices. Even though there is no applied voltage difference, a flow of 
electrons can tunnel from one superconductor to the other. This is due to the quantum 
mechanical waves in the two superconductors of the Josephson junction overlapping 
with each other. If we denote the difference in phases of the wave functions by (j) 
and the spatial and temporal variable along the junction by x and t, respectively, the 
electron flow tunneling across the barrier, i.e. the Josephson current, / is proportional 
to the sine of (j){x,t), i.e. / ~ sin(/)(x,i). In a long Josephson junction, the phase 
difference (j) satisfies a sine-Gordon equation. 
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In a standard long Josephson junction, the energetic ground state of the system 
is (pix) constant satisfying sin0 — 7, where 7 is an apphed constant (dc) bias current, 
which is taken to be zero here. A novel type of Josephson junction was proposed 
by Bulaevskii et al. [DIl], in which a non-trivial ground-state can be realized, char- 
acterized by the spontaneous generation of a fractional fluxon, i.e. a vortex carrying 
a fraction of magnetic flux quantum. This remarkable property can be invoked by 
intrinsically building piecewise constant phase-shifts 0{x) into the junction. Due to 
the phase-shift, the supercurrent relation then becomes / ~ sm{(f> + 9). Presently, 
one can impose a phase-shift in a long Josephson junction using several methods (see, 
e.g., and references therein). 

Due to this nontrivial properties, Josephson junctions with phase shifts may have 
promising applications in information storage and information processing [51 |H] . Be- 
cause of their potential applications, the next natural question is: what is the eigen- 
value of the ground state? It is important because Josephson junction based devices 
should not operate at frequencies close to the eigenfrequency of the system, as un- 
wanted parasitic resonances can be induced. 

The eigenfrequency of the ground state in the simplest case of Josephson junctions 
with one and two phase-shifts has been theoretically calculated in [71 [H [U [TOl HH [12] . 
More importantly, the eigenfrequency calculation in the former case has been con- 
firmed experimentally recently in [13l [14] . The experimental measurements were per- 
formed by applying microwave radiation of fixed frequency and power to the Josephson 
junction. At some frequency, the junction interestingly switches to the resistive state, 
characterized by a non-zero junction voltage. In terms of the phase-difference (j), the 
averaged Josephson voltage < > is proportional to 



where V is the domain of the problem and T ^ 1. It was conjectured that the driving 
frequency at which switching occurs is the same as the eigenfrequency of the ground 
state [13]. It is assumed that the jump to the resistive state is due to the resonant 
excitation of the breathing mode of the ground state and the applied microwaves, 
similar to the resonance phenomena observed in a periodically driven short (point- 
like) Josephson junction reported in [l5l [T6l [T7] . 

It was also noticed in [T3] that the accuracy of the microwave spectroscopy de- 
pends on the magnitude of the eigenfrequency. To measure a large natural frequency, 
the method requires an applied microwave with high power, which influences the mea- 
surement due to the nonlinearity of the system. Here, we consider an infinitely long 
Josephson junction with phase-shifts with no applied constant (dc) bias current. We 
show that in such a system, the breathing mode cannot be excited to switch the junc- 
tion into a resistive state, provided that the microwave amplitude is small enough. It 
is the case even when the applied drive frequency is the same as the eigenfrequency, 
because of higher harmonic excitations in the form of continuous wave emission. In 
other words, the breathing mode experiences radiative damping. Such damping is 
not present in short junctions as the phase difference (j) in that limit is effectively 
independent of x. This then confirms the observed experimental results. 

The governing equation we consider herein is given by 

(1.2) (t)^^{x,t) - (j)tt{x,t) ^ 8in{(j) + e{x)) + hcos{nt), x G M, t > 0, 

describing an infinitely long Josephson junction with phase-shifts driven by a mi- 
crowave field. Equation (jl.2l) is dimensionless, x and t are normalized by the Joseph- 



(1.1) 
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son penetration length Aj and the inverse plasma frequency lu~^ , respectively. The 
applied time periodic (ac) drive in the governing equation above has amplitude h, 
which is proportional to the applied microwave power, and frequency fl. Here we 
study two cases of the internal phase-shift 



(1.3) 9{x) 
with a < 7r/4, and 

(1.4) eix) : 



0, |a;| > a, 
TT, |a;| < a, 



0, a; < 0, 
— K, a; > 0, 



with < K < 27r, which are called — tt — and — k Josephson junction, respectively. 
These are the simplest configurations admitting a uniform and a nonuniform ground 
state, respectively. The phase field is then naturally subject to the continuity 
conditions at the position of the jump in the Josephson phase (the discontinuity), i.e. 

0(±a-) = </>(±a+), M±a~) = M±a+), 

for the — TT — junction and 

0(0-) = 0(0+), ^,(0-)- 0,(0+), 

for the — K junction. 

The unperturbed — tt — junction, i.e. (II. 2p and (|1.3p with h = 0, has 

$0 = 

(mod 27r) as the ground state and by linearizing around the uniform solution we find 
a breathing localized mode [2 



;os(aVl +w2)e^i~"'(''+^), x < -a, 
(1.5) $i(x,t) = e"^* { cos{xVru^), \x\ < a, 

cos(aVl + w2)e^i-'^'('^-^), x > a, 

with the oscillation frequency w given by the implicit relation 



(1.6) g = , ^ =tan~^ [ \ - — uj^ < 1. 

As for the unperturbed — k junction, i.e. (|1.2p and p.4p with h = 0, the ground 
state of the system is (mod 27r) 

^ / f 4tan-ie^«+^, a; < 0, 

<fo(a;,t) = | ^_4tan-ie---, a: > 0, 

where xq = lntan(K/8). Physically, $o (|l-7p represents a fractional fluxon that 
is spontaneously generated at the discontinuity. A scanning microscopy image of 
fractional fluxons can be seen in, e.g., [181 119] . Linearizing around the ground state 
(II. 7p we obtain the breathing mode [51 [TO] 

(1.8) $i(a:, t) - e | ^^(,^_,) ^^^^^^^^^ - a;) - A] , a; > 0, 
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with 



(1.9) 



and the oscUlation frequency 



(1.10) 



w(k) = ±W-cos- 




- + ^4-3cos^-j. 



In addition to the eigenfrequency (|1.6I) or (|1.10p . a ground state in a Josephson junc- 
tion also has a continuous spectrum in the range > 1. 

If a ground state is perturbed by its corresponding locahzed mode, then the 
perturbation will oscillate periodically . The typical evolution of the initial condition 



for some small initial amplitude Bq and h — is shown in the top and bottom panel 
of Fig. 11.11 for the two cases above. For a — tt — junction, one can see a clear 
mode oscillation on top of the uniform background state = 0. In the case of a 
— K junction, the periodic oscillation of the localized mode makes the fractional kink 
oscillate about the point of discontinuity, a; = 0. 

Using a multiple scale expansion, we show that in the absence of an ac drive, such 
a breathing mode oscillation decays with a rate of at most 0(t~^/'^) and 0{t~^/'^) for 
a junction with a uniform and nonuniform ground state, respectively. The coupling 
of a spatially localized breathing mode to radiation modes via a nonlinearity and the 
same decay rates have been discussed and obtained by others in several contexts (see 
PP] and references therein). Interactions of a breathing mode and a topological kink, 
creating the so-called 'wobbling kink', or simply 'wobbler', have also been considered 
before, see [IHIISI] (see also references therein) for (j)^ wobblers and [351 [521 Hi] for sine- 
Gordon wobblers. Nonetheless, the problem and results presented herein are novel 
and important from several points of view, which include the fact that our fractional 
wobbling kink is in principal different from the 'normal' wobbler. Usually, a wobbler 
is a periodically expanding and contracting kink, due to the interaction of the kink 
and its odd eigenmode. Because our system is not translationally invariant, our 
wobbler will be composed of a fractional kink and an even eigenmode, representing 
a topological excitation oscillating about the discontinuity point (see also [24] for 
a similar situation in discrete systems, where a lattice kink interacts with its even 
mode). Such an oscillation can certainly be induced by a time periodic direct driving, 
as considered herein. More importantly, our problem is relevant and can be readily 
confirmed experimentally (see also, e.g., [IHl 123 [2H] for experimental fabrications of 
— TT — Josephson junctions). 

The presentation of the paper is outlined as follows. In Section 2, we construct a 
perturbation expansion for the breathing mode to obtain equations for the slow-time 
evolution of the oscillation amplitude in a — tt — junction by eliminating secular 
terms from our expansion. In Section 3, the method of multiple scales is applied to 
obtain the amplitude oscillation in the presence of driving, extending the preceding 
section. In Sections 4 and 5, wc apply the perturbation method to the wobbling kink 
in a — K junction. We confirm our analytical results numerically in Section 6. We 
also show in the same section that there is a threshold drive amplitude above which 
the junction switches to the resistive state. Yet, we observe that the switching to the 
resistive state is due to the modulational instability of the background. We conclude 
the present work in Section 7. 



(1.11) 



(^ = $o(a;)+So$i(x,0) 
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Fig. 1.1. The typical dynamics of a breathing mode (top) and a wobbling kink (bottom) in an 
undriven — tt — and — k junction, respectively. 

2. Freely oscillating breathing mode in a O — TT — O junction. In this section 
we construct a breathing mode of tlie sine-Gordon equation (|1.2p with h = and 9 
given by (|1.3I) . 

We apply a perturbation method to the equation p.2p by writing 

(2.1) = 0o + e0i + e^'/'2+e^<?!'3 + ---, 

where e is a small parameter. We further use multiple scale expansions by introducing 
the slow time and space variables 

(2.2) X„ = e"a;, T„ = e"t, n = 0, 1, 2, . . . , 

which describes long times and distances. In the small limit of e, the different scales 
above become uncoupled and may be considered as independent variables. 
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(2.3) 



In the following, we use the notation 

d 



d„ 



such that the derivatives with respect to the original variables in terms of the scaled 
variables using the chain rule are given by 



(2.4) 
(2.5) 



d_ 

dx 
d_ 



Do + eDi + e^D2 + e^D3 



Substituting these expansions into the perturbed sine-Gordon equation (II. 2p along 
with the expansion of (j) and equating like powers of e, we obtain a hierarchy of linear 
partial differential equations; 



(2.6) 
(2.7) 



0(1) : d^(bo-Dl(bo=smi0 + M, 

0{e) : 5201 - Dl(f>i - cos(0 + 0o)0i = 2i?o^i0o - 2dodi^o. 
Solutions to the equations above for the — tt — junction are given by 

(2.8) 0o(^o,ro) = O, 
and 

(2.9) </)i(Xo,ro) = B(Xi,..,Ti,..)$i(Xo,ro) + c.c., 

where $1 is given by (jl.Sp . B{Xi, . . . , Ti, . . .) is the amplitude of the breathing mode, 
which is a function of the slow time and space variables only. Throughout the paper, 
c.c. stands for the complex conjugate of the immediately preceding term. 
Next, we consider the 0{e^) equation 

(2.10) d^(j)2 - Dl(j>2 - cos{e + 0o)02 = 2DoDiq)i - 2dodi(j>i. 
Evaluating the right hand side for the different regions, we obtain 

,Vl-w^(Q+^o)+iwTo 



doh - Dl(j)2 ~(j}2^ 2cos(aVl +^2) iujDiB - vT^Z/^diB 
dl(t>2 - Dl(t>2 + 02 2 \iujDiB cos{Xo^/T^^) + ^/l + uj'^diBsui{Xo^/l + uj"^) 



dl(j)2 - Dl(j>2 - 02 = 2cos(aVTTw2) \iujDiB + yJl-u^diB 



for Xq < —a, \Xq\ < a, and Xq > a, respectively. These are linear wave equations 
with forcing having frequency uj. 

Substituting the spectral ansatz (t)2{XQ,To) = 02(-'^o)e'"^°, we obtain the cor- 
responding set of ordinary differential equations with forcing term which has the 
frequency uj 



^002 - (1 - ^2)02 = 2 cos(aVT+^) iioDiB - yj\ - uj^diB 



5o02 + (1 + C^' 



iujDiB cos{XaVT+u^) + ^/T+u^^lB sm{XoVT+u^) 



dl4>2 - (1 - t^^)02 = 2 cos(aVl+^) iujDiB + \/l - u'^diB 
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We need to find the bounded solution of ttie above equations, wliicli are of tlie form 
Tip (Xq) = f (Xq) where T is a self-adjoint operator given by the left hand side of the 
above system. 

The Fredholm theorem states that the necessary and sufficient condition for the 
above nonhomogeneous equation to have a bounded solution is that its right-hand 
side /{Xq) be orthogonal to the complete system of linearly independent solutions of 
the corresponding homogeneous equation, Tip(Xo) = 0. By applying the theorem, we 
find the solvability condition 

(2.11) DiB = 0. 

The bounded solution of (|2.10p is then given by 



62 = diBe 



ibjTo 



^^^^X q (j^^ _ COS {aVTTuP) e ^/i^(° +^o) -|- c.c, Xq < -a, 
(XqVTT^) C22 - Xo COS (XpVl + ^2) ^ c^c^^ ^Xo\ < a, 

-VT^Xo Q^^^ _ x„ cos {aVT+uP) evT^("-^o) + ^.c, Xq > a, 



e 

cos 

e 



where C21 = C'23 and C22 are constants of integration that have to be found by 
applying the continuity conditions at the discontinuity points X^ = ±a. 

It should be noted that diB as well as 9„i3 in later calculations do not appear in 
the solvability condition. Therefore, we take the simplest choice by setting 

diB = 0. 

This choice is also in accordance with the fact that if diB was non zero, then {e^(t>2) 
would become greater than (e^i) as Xq — >■ ±00 due to the term [XQe^^~^'^^'^'^'^°^) in 
the expression of 02 above, leading to a nonuniformity in the perturbation expansion 
oicf). 

The equation at the third order in the perturbation expansion is 

(2.12^2^3 - Dlcp^ - cos(0)(/)3 = 2{DoD2 - 8082)^1 + {D} - 8f)cPi - ^q^l cos(0). 



Having evaluated the right hand side using the functions (j)o and 0i , and splitting the 
solution into components proportional to simple harmonics we get 



Jon 



Fi, Xo< -a, 
F2, \Xo\<a, 
F3, Xo > a, 



where Fi, F2, F3 are given by 



Fi = 



2 (tLjD2B - Vl - w^aaS) cos (a^l+w^) e^^^'^''+^"^+''^^° 
- -B\B fcos'' (aVT+^) e3^i=^('^+^»)+^-^« 



-is^ cos3(av/m^)e3^i^(''+^«)+3*-^«. 



Fo - 



2 iijjD2B cos ( y/Y+^Xo ] + 282B Vl + u)'^ sin 



islSpcos^ (^l+uj-^Xo 



-B^ cos- 
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F3 = 2 (^iujD2B + yj\ - Lj^d2B^ COS (aVTT^) e^''i^('^--^o)+*'^^« 
^Ib\B\^ cos' (a^/lTI^) e3^T^(«-A-o)+.-To 
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These are linear wave equations with forcing at frequencies lu and 3a;. The former 
frequency is resonant with the discrete eigenmode and the latter is assumed to lie in 
the continuous spectrum (phonon band), i.e. 

(2.13) 9uj^ > 1. 

As (I2.12p is linear, the solution can be written as a combination of solutions with 
frequencies present in the forcing, i.e. 



(2.14) 



^3 + ^3 



+ c.c. + 



,(2) 2iwTo 



+ C.C. + 



A3) 3iuiT„ 



+ C.C. 



This implies that ^3 satisfies the following nonhomogeneous equations in the three 
regions below 



2iujD2B cos {aVTTo?) eVi-""('^+^o) 

B\B\^ cos3 ( gyTT^ ) e' ^^i-+Xo) ^ Xo < -a, 

ao^</>« - (cos(0) - = 2 ^^D^Bcos if+^o) 

°' V ^"^3 ) +lB\B\^ cos' (Vl + w^Xo) , |^o|<a, 

2iujD2B cos (aVl + w2) f.VT^(a-x„) 
-\B\B\^ cos' (aVl +a;2) e3^^i-"'(''--^")^ > a. 
In the above equation, it should be noted that we have imposed the fact that 

d2B = 0, 

as we previously discussed. 

The solvability condition for the first harmonic gives 

D2B = ikiB\B\^, 



(2.15) 

where 

(2.16) 



ki = 



(3 - 7 _ 2 _ 2 a;2 + 6 Vr 



cj* tan 



32 (l + cj2 + VT^7J4tan-i (y^i^)) 



The solution for the first harmonic is then given by 



(2.17) 

where 



vi{Xo), Xo < -a, 
'^'^\Xo,To) = B\Bf { V2iXo), \Xo\<a, 
vsiXo), Xo > a. 



viiXo) = c^^''°C3i - 



v/2(l+L^2) |^fc3^g2 vT=^Xo _ kl^^ ^-VT^Xo 



64 w y/T~~uj^ Ui 



V2{Xo) = Re{^ 



/c32Re(^ 



- fc32lm(e 



iVi+<^^Xo 



V3{Xo) = e'^^''«Q 



33 



16 VI +W2U2 

32 V(2-a;2)(l 
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The expression of the functions fcsj and /csj, j = 1,2,3, is given in (jA.l[) - (jA.6p . The 
coefhcients C31 — C32 and C33 are constants of integration, that should be determined 
from the continuity conditions at = ±a. 

We do not consider the equation for the second harmonic <i>\ as it does not 
appear in the leading order of the sought asymptotic equation describing the behavior 
of breathing mode amplitude. 

The equation for the third harmonic (/ig is 



[ -^cos^( N/rT^ a)e3^'^^^(°+-^°\ Xo<-a, 
dl<i^ - (cos(0) - 9c^2) 0(3) ^ J 1 cos^V T+Z^X q), \Xo\ < a, 

[ -icos3(Vl +w2a)e3^''T^('^-^«), Xq > a, 

whose solution, using the same procedure as above, is given by 

r P3i(^o), Xo<-a, 
0f (Xo,To) = b3<^ ^32(^0), |^o|<a, 
[ P33(^o), ^0 > a, 

where 

P32(Xo) = cos{Vl + 9^Xo) Caa - (w' - 3) cos (Xo VTT^) , 

P33(^o) = e-^^^^°C;3 - i cos3 (aVl+^) e^^^(-^o)^ 



and C31, C32, and C33 are non-zero constants of integration that also are determined 
from the continuity conditions at the discontinuity points. 

Note that due to the assumption (|2.13p . the second term in P3i{Xq) and P33(Xo) 
will decay to zero while the first term tends to C31 and C33 as Xq — > Too, respectively. 
This implies that 6'^'''^'^° (f)'^'' +c.c. represents a continuous wave radiation traveling to 
the left and right. 

Solving the O (e^) equation we obtain 

(1)4 = 0, 

from the solvability condition 

(2.18) D3B = 0, 

which is similar to the case of 4>2 ■ 

Equating terms at C(e^) gives the equation 

5o'<^5 - ^o'-^s - cos(0)05 = 2iDoDi - dodi)^i + 2{D^Di - d^di)4>i + (Dj ~ dl)4>i 

-,2 a2\ - , o/n n a » ^^ , / ^ j, 2 , , 1,5 
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Having calculated the right hand side using the known functions, we again split the 
solution into components proportional to simple harmonics as we did before. The 
equation for the first harmonic is given by 



(2.19) 



^2^,(1) 



Gi, Xq < —a, 



{cosi0) - Lu^) = { G2, \Xo\<a, 
G3, Xq > a, 



+ e 



where 

Gi = 2iuj D^B cos (^a\/l 

-B\B\^ fci^cos (av/Tw) e^^^^''+^"'> + 2ujkiVi{Xo) 
+ i cos^ (av/l+c.2) e^^^^(-+^»H3v,iXo) + P3iiXo)) 

G2 =2ia;L>4Bcos(^XoyTT^) - B\B\* fci^cos (Xo\/l + uj^^ + 2u}kiV2{Xo) 
-i cos^ (3 1^2(^0) + ^32(^0)) + ^ cos^ (Xo VTT^) 

G3=2iuj D4B cos (a VTTI^) e^^^^''-^'>'> 

-B\B\^ [/ci^ cos (av/TT^) e^''i^(''-^") + 2cjfcit.3(Xo) 



cos 



2] „2Vl^(a-Xo) 



(3u3(Xo) + i'33(^o)) 
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COS 



Here, vi{Xq), V2{Xq), V3 (Xq) are the bounded solutions of (I)''^\Xq, Tq) and P3i{Xo), 

P32{Xo), ^33(^0) are the bounded solutions of (j)''r^\Xo,To) as solved above. 
The solvability condition of (I2.19p is 



(2.20) 
where 



T2^ 
'2uj'^{uj] 



D4B = k2B\B\\ 



e C, T2 = fc?*(t^) + 2ujkiaXo) + a{Xo) + piXa) + 7(^0), 



*(w) 
C(^o) 



(x/TT7^ + ^/T^tan-i (/iiir)) 



t;i(Xo) cos (av/l+w2) e^^'i^("+^'')dXo + j ^2(^0) cos (XoVTT^) (i^o 
1^3(^0) cos [a^/T+^j e^i^("-^«)dXo, 
^(^o) = i/ j3«i(Xo) + P3i(^o))cos3 (av/rT^)e^^(''+^"'dXo- '"'^^^^^^ 
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fi{Xo) (3^^2(^o) + ^^32(^0)) cos^ (Xa^/TV^) dXo + ^ j co&\Xo^/T+J^)dXo, 



72V1 - a;2 



We shall postpone the continuation of the perturbation expansion to higher orders 
as we have obtained the oscillating and decaying behavior of the breathing amplitude 
(|2.15l) and (I2.20p . which is our main objective. 

By noting that 

dB 

— = eDiB + e^D2B + e^DgS + e^D^B, 
dt 

and defining b = eB, ([2lT|) . (|2?T5l) . (|2?T8l) . and (fl!20|) can be combined to obtain 

(2.21) — = kMbf I + k2b\b\\ 

dt 

One can calculate that the solution of (|2.2ip satisfies the relation 



(2.22) mi 



1- 4:Re{k2)\b{0)\H 



where 6(0) is the initial amplitude of oscillation. By assuming that Re(fc2) < 0, as 
will be shown later, this equation describes the gradual decrease in the amplitude of 
the breathing mode with order ©(t^^/^) as it emits energy in the form of radiation. 

Remark 2.1. The O decay of the oscillation amplitude is because of our 

assumption 112.13]) . If one has (Scj)^ < 1 instead, then the decay rate will be smaller 
than O (t^^/^^ as the coefficient k2 in 112. 20\) will be purely imaginary. This then leads 
us to the following conjecture. 

Proposition 2.2. If n > 3 is an odd integer such that 

l/{n-2f >uj^ > l/n^, 

then the decay rate of the breathing mode oscillation m — vr — Josephson junctions 
with a < 7r/4 is of order 0(t-i/(«+i)). 

This proposition implies that the closer the eigenfrequency uj to zero, ie. a — )■ 7r/4, 
the longer the life time of the breathing mode oscillation. 

3. Driven breathing mode in — tt — junction. We now consider breathing 
mode oscillations in a — tt — junction in the presence of external driving with 
frequency near the natural breathing frequency of the mode, i.e. (|1.2p and (|1.3p with 
h^O and = w(l + p). 

By rescaling the time fit — ujt, (|1.2I) becomes 



(3.1) 0,,(x, r) - (1 + p)^cj)rr{x, t) = siu {4> + 9) + hi (e'"" + c.c.) . 

Here, we assume that the driving amplitude and frequency are small, namely 

(3.2) h = e'H, p = e^R, 
with H,R^O{l). 
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Due to the time rescaling above, our slow temporal variables are now defined as 
(3.3) X„ = e"a;, r„ = e"r, n = 0,1,2,.... 



In the following, we still use the short hand notations 

Performing a perturbation expansion order by order as before, one obtains the 
same perturbation expansion up to and including order as in the undriven case 
above. At third order, we obtain 

(3.4) 92^3 - Dl(l)z - cos(0)(/.3 = {Dl - )(/.! + 2{DoD2 - dod2)(bi - ^(/-f cos(0) 

+ (e-- + c.c.) . 

The only difference from the undriven case is the presence of a harmonic drive in the 
last term. 

The first harmonic component of the above equation gives us the solvability con- 
dition 

(3.5) D2B = kiB\B\^ i + hHi, 
where 



(3.6) h = , 

V2io [1 + c^2 + yi^T^tan-i (^t+^)) 

and ki is given in (j2.16p . The solution for the first harmonic can then be readily 
obtained as 

, B\B\^vi{Xo)+Hv\{Xo), Xo<-a, 

) u\ nl2 



where 



^~i(Xo) = e^^^"H3i 



B\B\'v2iXo) + Hv2iX„), \X„\ < a, 
B|B|2w3(^o) + HvsiXo), Xo > a. 



e 



2 (l-w2)3/2^^ 
u+Xo ^(^4 _ ^-^Xo + VT~~^{1 + UJ^)) - 2 \/r^(l + ^2) 



V2iXo) = Rc(e*^''i+^^'')H 



4 

1 



32 



2(1+^2) 

Vl + a;2XoIm (e* ^^^^-^o) + Re(e^ VT+^^o) 



V2VT 



up- U2 



v,{Xo) = e-^^"H33 

2 ui 



e-^'^i^-^»((e" - 2e^^'i^^«)\/l-a;2(i + a;2) _^ 2 e''Xo(l - w^)) 

4(l-w2)3/2y^ ' 

and tt = \/r~~cj2 tan~^ ( ^+1;""^ ) • '^^^ constants of integration H31 = H33 and H32 
are determined by applying the continuity conditions at the discontinuity points. 



Breathing modes of long Josephson junctions with phase-shifts 13 

One can check that the solution for the third harmonic ^3'^^ (Xq, Tq) is the same as 
in the undriven case as well as ^3 , which has no role in the leading order asymptotic 
expansion. 

The equation at 0{e'^) is 

5o^04 - £'o'<^4 - cos {0 + 0o) 04 = 2 {DqDi - dodi) cf>3 + 2 {D1D2 + DoDs) 0i 

-2 {did2 + dod3) (j>i +4RDo^(j)i, 

with the solvabihty condition for the above equation 

D3B = -2iujBR. 

This impUes that ^4 = as for the case of 02- 
At C(e^), we obtain 

SoVs - - (As cos(e) = 2{DoDi - dodi)(l>i + 2{DiD^ - d3di)(f>i + {D^ - dl)(f>i 

+{D! - a?)<A3 + 2(£>2-Do - d2do)h - ^i^h - Y^<^i') cos(^) +iRDoDi(f>i. 

Evaluating the right hand side, we again spht the solution into components propor- 
tional to simple harmonics as we did before. For the first harmonic, we obtain that 

r Ml, Xo < -a, 

(3.7) SoVi'^ - (cos(e) - = <^ M2, \Xo\ < a, 

[ Ms, Xo > a, 

where 

Ml = (2iujD^B- ki^B\B\'^ - kxh\B\^H^ cos (a^/l+^^ ^VT^{a+Xo) 

-2 w {kiB\Bt + h\B\'H^ vi{Xo) + B\B\'' cos^ (a\/T+^) v'^^z^C^+^o) 

B\Bf (3t;i (Xo) + P31 (Xo)) cos' (aVl+w') e'-^I^(«+^'') 
-\h{2 \B\^ + B^) (Xo) cos' (a^/^+^) e^^/^(-+^'>) ^ 

M2 = {2iLuD^B - ki^B\B\^ - kih\B\^H) cos (XoVTTo^) 

-2uj {kiB\B\^ + h\B\''H)v2{Xo) -^B\B\^cos^ (Xo ^1 + w') 
+^ (3t;2 (Xo) + P32 (Xo)) cos' (XoVl + w^) 

+ i if (2 + B') ^72(Xo) cos' (XoVl+^') , 

Ms = (2ioj D4B - ki^B\Bf - kihlB^H^ cos (ay/l + w') e^'^^^""^") 

-2 w (fciB|B|* + ZilS]'//) W3(Xo) + ^ cos^ (a Vl + e^^-^i^^""^") 

-ii?|i3|^(3t;3(Xo) + P33(^o))cos' (ayiT^) e^v^Ca-Xo) 
-1 H (2 \Bf + B') {}3(Xo) cos' (a^TT^) e' ^^^(-^0). 



14 



ALI, SUSANTO, WATTIS 



The solvability condition for (|3.7p is 

(3.8) D4B = fc2 B\B\'^ + Hi {l2\B\^ + hB^) 

where li is given by p.6p and 

'2 = ^ .T. / M '3 



2a;*(a;)' 2tj^'(w)' 
= fci?i*(c.) + 2t^?iC(^o) + 2 («2(^o) + /32(^o) + 72(^0)) , 
T3^2 - a2(Xo) + /32(^o) + 72(^0), 

«2(^o) = ^ I %~i(Xo)cos3(aVT+^)e3^^(°+''"^d^o, 

/32(Xo) - / ^2{Xo) cos^ (Xo VTT^) dXo, 

72 (^0 ) = 2 y ) (« v/l+^) ^/^('^-^"^ dXo . 

So far we have obtained the sought leading order behavior of the breathing amplitude. 
Performing the same calculation as in p.2ip . we obtain the governing dynamics of 
the oscihation amphtude in the presence of an external drive 

(3.9) '^'^^^kih\b\'^i + hhi + k2b\b\^ + h{l2\b\'^ + hb'^) i-2iu;bp. 

One can deduce that a nonzero external driving can induce a breathing dynamic. It 
is expected that for large t, there will be a balance between the external drive and 
the radiation damping. 

4. Freely oscillating breathing mode in — k junction. In this section, we 
will consider (|1.2p with 9 given by (|1.4p . describing the dynamics of the Josephson- 
phase in the — k long Josephson junction. 

By applying the method of multiple scales and the perturbation expansion as 
before to the governing equation, we obtain from the leading order 0(1) and 0{e) 
that 

(4.1) ^o = MXo), 01 -s(Xi..,ri..)$i(Xo,ro) + c.c., 

where $0 and (f>i are given by (II. 7p and (II. 8p . respectively. 

Using the fact that 0o is a function of Xq only, the equation at O(e^) is 

(4.2) 92^2 - Dlcj)2 - cos{e + = 2DoDicj)i ~ 29o9i</'i - y sin((? + 

After a simple algebraic calculation, one can recognize that the right hand side of 
the above equation consists of functions having frequencies 0, w and 2lu. Therefore, 
solutions to the above equation can be written as 

(4.3) 02 = 4°^ + 4'^e^"^» + c.c. + 0f e2'"^« + c.c, 

which implies that (j)^"^ , and 02^'' satisfy the following nonhomogeneous equations 
in the two regions Xq < and Xq > 



'9o2-cos((? + 0o))0f 



e 



tanh(a;o + Xq) — \/l — cj^ 



tanh(a;o — Xq) — yl 



2 

sin(0o), -'^o < 0, 

2 

sin(0o ~ 1^); Xq > 0, 
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jDiB - diB^/l-Lo^) (tanh(a;o + Xq) - ^JI~lo^) 



— 9ii?sech [xq + Xf)) 
2 



Xo < 0, 



[iLoDiB + diBVl - uj'^) (tanh(a;o - Xq) - 
-9ii?sech^(xo — Xq) 



Xn > 0, 



' 4a;^ - cos 



(2) 



B^ 
2 



tanh(a;o + -^o) — \/l — 



tanh(a;o — Xq) — 



sin(0oj, 



^0 <0, 



sin(0o — k), Xq > 0. 



By using arguments as in the preceding sections, we set 

DiB = 0, diB = 0. 

Hence, we find that 

4'\Xo,To)=0. 
The solution for the other harmonics are 



(4.4) 
(4.5) 



4°) = 

a(2) 



B\ 



= B' 



Eo{Xo), Xo<0, 
io(Xo), Xo>0, 

E2{Xo), Xo<0, 
MXo), Xo>0, 



where Eo{Xo), Eq{Xq) ,E2{X^) and ^2(^0) are given in (lA?7l)- (lAlO| . Cqi, C02, C21 
and C22 are constants of integration that should be found by applying continuity 
condition at the point of discontinuity Xq — 0. 

In the following, we will assume that the harmonic 2a; is in the continuous spec- 
trum, i.e. 



(4.6) 



Auj^ > 1. 



With this assumption, (j)^\XQ,TQ) will not decay in space and e^*"^''(/)2^'' (Xq, Tq) 
describes right moving radiation for positive Xq and left moving radiation for negative 

Equating the terms at O(e^) gives the equation 



5^(/)3 - 1^0-^3 - cos( 



h)h = (2DqD2 - 2dQd2)<Pi + (Di ~ dt)^i + 
{2DqDi - 2dQdi)4>2 - 4>ih sin(0 + (Pq) 



6 

where we have used the fact that (j>Q depends only on Xq . Having calculated the right 
hand side using the known functions 0o, 4>i and 02, we again split the solution into 
components proportional to the harmonics of the right hand side. Specifically for the 
first harmonic, we have 



(5o^+w2_cos(6l + 0o))4 



(1) _ / Li, Xq < 0, 
L2, Xq>0, 
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where 

Li = 2iujD2B 



t anil (.To + Xq) — \/\ — up- 



,VT=i;j^(a;o+Xo) 



~B\Bf {Eo{Xo) + E2{Xo)) sin{^o) 
-^B\Bf cos(?!>o) tanh(a;o + Xq) - y^l-co^ 



tanh(.To + Xa) - VT 



„3 VT^{xo+Xo) 



L2 = 2ioj D2B 



tanh(xo — Xq) — \/l — u)'^ 



~B\B\''(e^{Xo)+E2{Xo, 



sin( 00 ^ 



tanh(a:;o — Xq) — \/ \ — ixP- 



— \^\B^ cos(0o — tanh(xo — Xq) — \/\ — uP- 
The solvabihty condition of the equation will give us 

where 

Ti 

mi = — 



,3 ^/T^{xo-Xq) 



/i(^o) 



hiXo)dXo+ f2{Xo)dXo + - h{Xo)dXo + - f4{Xo)dXo, 

-00 Jo ^ J -00 ^ Jo 

2oje'^^'^^^° (2 VI - + 2 - _ e2^o (2^/1 - + - 2) ) 



\/l-w^(l + e2^o) 
-2 (£;o(^o) + £^2(^0)) sech (xq + Xq) tanh(xo + Xq) 

[tanh(:Eo +^0) - Vl^] e'^^^^'''+^"\ 
f2{Xo) = 2 (^o(^o) + ^2(^0)) sech (xo - Xq) tanh(xo - Xq) 



tanh(a;o — Xq) — \/\ — iJ^ 



JlsfT^{xo-Xo) 



fs{Xo) 
h{Xo) 



tanh(a;o +Xo) - ^1-^2 (1 - 2sech2(xo + Xo))e^^^^^(^''+^°\ 



tanh(a;o - Xq) - Vl - (1 - 2sech2(xo - Xo))e^^ 



We will not continue the perturbation expansion to higher orders as we have obtained 
the leading order behavior of the wobbling amplitude. 
Using the chain-rule and writing b = eB, we obtain 



(4.7) 

It can be derived that 



at 



2Re(TOi)|6|^ 



with the solution given by 
(4.8) |6| = 



\b{OW 



l-2Re(TOi)|6(0)|2t' 
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and the initial amplitude 6(0). It can be clearly seen that the oscillation amplitude of 
the breathing mode decreases in time with order 0{t-'^^^). 

Remark 4.1. Similarly to Remark \2.1{ the 0{t^^^^) amplitude decay is caused 
by the assumption I14.6^ . One therefore can introduce a similar proposition as before. 
This then leads us to the following conjecture. 

Proposition 4.2. If n> 2 is an integer such that 

l/{n-lf >L0^ > l/n^, 

then the decay rate of the breathing mode oscillation in — k Josephson junctions is 
of order 0{t~^/"). 

5. Driven breathing modes in a — k junction. We now consider breathing 
mode oscillations in — k junctions in the presence of external driving with frequency 
near the natural breathing frequency of the mode, i.e. (|1.2I) and (|1.4p with h and 
ft = + p). Taking the same scaling as in the case of driven — tt — junction, 
we obtain (13.11) . Yet, here we assume that the driving amplitude and frequency are 
small, i.e. 

(5.1) h = e^H, p = e^R, 
with H,R^O{l). 

Performing the same perturbation expansion as before, up to 0{e) we obtain the 
same equations as in the undriven case, which we omit for brevity. 
The equation at O(e^) in the perturbation expansion is 

(5.2) 90^2 - Dl^2 ~ cos{d + (bo)cf>2 = 2 {DoDi - dodi) 0i " y sm{0 + ^ 

+ ii/(e-- + c.c.). 

Again, one can write the solution (j32 as a combination of solutions with harmonics 
present in the right hand side. In this case, the first harmonic component is different 
from the undriven case due to the driving, yielding the solvability condition 

(5.3) DiB = mHi, 
where 

r]{xo,uj) 



m 



dXo 



T]{xo,uj) = / e^''^^(™+-^'') tanh(2;o + Xq) - ^1-^2 

J —OO 

+ / e^^i^(""-^«) tanh {xg - Xq) - dXo. 
Jo 

The solution for the first harmonic is then 

(5.4) 0, (Xo,To)-i/| ^^^^^^^^^^^^^ ^^^Q^ 
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where 

(5.5) g(Xo) 

(5.6) ~g{Xo) = 

(5.7) n{Xn) 



, vixo,Xo) 

(2 V]~7J2 - 2 + cj2 + a;2e-2 i^o-Xo)-^ ^VT^(xo-x„)(j^^ 



(2 + 2-UJ^- Cj2e2 (^o+Xo)) g^/T^(.o+Xo)c;l 

I _|_ g2(2;o+^o) 

(^2g2 (xo+Xo) _ 2 yiT- 7J2 _ 2 + a;2) e^''i^(^o+^«Ui(Xo) 



(5.8) n(Xo) 



(^2g2 ix„+Xo) + 2 VT^^ - 2 + cj2-) e-vT^(xo+Xo)^2(Xo) 
4a;Vl-c^2(i + e2(^o+Xo)) ^ 
(^2g-2 (xo-Xo) ^ 2 yr^7J2 - 2 + cj2^ e^i^(^"-^«))a2 



(^2g-2(^„-Xo) _ 2 ViTr^ + 2 - a;2) e-^'^^(^«-^"M3(Xo) 

4w4Vl - ^2(1 + e-2(^o-Xo)) 
(g-2 (a;o-Xo)^2 _^ 2 VI - u;2 - 2 + a;2)e^''i^(^"--^«) A4(Xo) 
4a;4VT^^(l + e-2(^o-Xo)) 

Here, 51 (Xq), 51 (-'^o), and Aj, j = 1, . . . ,4, are given in (|A.lip -(|A.16p. Cgi, Cg2, Chi 
and Ch2 are constant of integration chosen to satisfy continuity conditions. The other 
harmonics are the same as in the undamped, undriven case. 
Equating the terms at 0{£^) gives the equation 

duct's - Dl4>3 - cos(0 + (/)o)03 = 2{DoD2 - doO^)^! + {Dj - df)(t>i + 2(DoDi - dodi)4>2 

+2RDl(j)i - (j)i(t)2 am{e + 0o) - y cos(6' + 0o)- 

The solvabihty condition for the first harmonic of the above equation gives 

(5.9) D2B = miB\B\^ -uBRi. 

Equations (|5.3I) and (15.91) are the leading order equations governing the oscillation 
amplitude of the breathing mode. Combining both equations gives us 

(5.10) — — = mhi + mi6|6P — iujbp. 

Q at 



Similarly to p.9p . one can also expect to obtain the result that a nonzero external 
drive amplitude induces breathing mode oscillation. 

6. Numerical calculations. To check our analytical results obtained in the 
above sections, we have numerically solved the governing equation (jl.2p with 9[x) 
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Fig. 6.1. Oscillation amplitude of the breathing mode in a — it — (top panel) and — re 
(bottom panel) junction. The solid curves are from the original governing equation ( clearly 
indicating the decay of the oscillation. Analytical approximations \2.22'i for the top panel and \4.8'^ 
for the bottom panel are shown in dashed lines (see the text). 



given by (jl.Sp or (ll.4p . We discretize the Laplacian operator using a central differ- 
ence and integrate the resulting system of differential equations using a fourth-order 
Runge-Kutta method, with a spatial and temporal discretization Aa; = 0.02 and 
At = 0.004, respectively. The computational domain is x G {—L,L), with L = 50. 
At the boundaries, we use a periodic boundary condition. To model an infinitely 
long junction, we apply an increasing damping at the boundaries to reduce reflected 
continuous wave incoming from the boundaries. In all the results presented herein, 
we use the damping coefficient 

, . ( {\x\ -L + Xa)) /Xa, \x\ > [L - Xa), 

^ ' " - \ 0, |x| < (L-x„), 

i.e. a increases linearly from Q; = Oata:: = ±(L — a;Q,)toa = lata: = ±L. We have 
taken Xa =20. To ensure that the numerical results are not influenced by the choice 
of the parameter values above, we have taken different values as well as different 
boundary conditions and damping, where quantitatively we obtained relatively the 
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same results. 

In this section, for the — tt — junction we fix the facet length a = 0.4, which 
implies that uj « 0.73825, and for the — k junction we set k = n, which implies that 
xq « —0.8814 and w ss 0.8995. For the choice of parameters above, we obtain the 
coefficients in the analytically obtained approximations (I2.2ip . (13. 9p . (|4.7p . and (|5.10p 

as 

0.0433, fc2 = -0.00324 -0.01400 i 

(6.2) ^1 = 0.6068, /2 = -0.10027, = 0.05934, 
TO = -0.6236, mi = -0.0182 - 0.0809 i. 

First, we consider the undriven case, h — 0. With the initial condition (jl.lip and 
Bq particularly taken to be 

(6.3) Bo = 



*i(0,0)' 



where ^i{x, t) is given by ()1.5p for — vr — junctions and (|1.8p for — k junctions, we 
record the envelope of the oscillation amplitude (/'(0,t) from the governing equation 
(|1.2p . In Fig. 16. 1[ we plot in solid lines (f){Q. t) of the — tt — and the — k junction 
in the top and bottom panel, respectively. 

From Fig. 16.11 one can see that the oscillation amplitude decreases in time. The 
mode experiences damping. The damping is intrinsically present because the breath- 
ing mode emits radiation due to higher harmonics excitations with frequency in the 
dispersion relation. 

It is then instructive to compare the numerical results with our analytical calcu- 
lations. With the initial condition 

6(0) = Bo/F, 

the analytical approximation is then given by F\h{t)\, where \b{t)\ is given by (|2.22p 
and (|4.8p for the — vr — and — k Josephson junctions, respectively. 

In general, the factor F is simply F = 2$i(0,0), i.e. F = 2 for the — tt — 
case and F k, 1.56 for the — k case. Yet, by treating as a fitting parameter we 
observed that the best fit is not given by the aforementioned values. For the initial 
condition (j6.3l) . we found that an optimum fit is respectively provided by = 2.06 and 
F — 1.8. Shown in dashed line in Fig. 16.11 is our approximation, where one can see a 
good agreement between the numerically obtained oscillation and its approximations. 
In the top panel, the approximation coincides with the numerical result. 

Next, we consider the case of driven Josephson junctions, i.e. (II. 2p with h ^ 0. 
In this case, the initial condition to the governing equation (II. 2p is (|l.lip with Bq 
particularly chosen to be 

Taking specifically ft = uj, respectively we present the amplitude of the oscillatory 
mode (f>{0, t) of — tt — junctions with h — 0.002 and h — 0.003 in the top and 
middle panel of Fig. 16.21 These are the typical dynamics of the oscillation amplitude 
of the breathing mode, where for the first case the envelope oscillates periodically in 
a long time scale and for the second case the amplitude tends to a constant. 

To compare it with the asymptotic analysis, we have solved the amplitude equa- 
tions p.9p and (I5.10p numerically using a fourth-order Runge-Kutta method with 




Fig. 6.2. The same as in Fig. [O but for nonzero driving amplitude. Top and middle panel 
correspond to driven — tt — junctions uiith h = 0.002 and h = 0.003, respectively. Bottom panel 
corresponds to driven — k junctions with h = 0.01. 
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a relatively fine time discretization parameter, as exact analytical solutions are not 
available. The analytical approximation is again given by where F in this 

case is taken to be exactly F = 2$i(0,0). It is important to note that ideally p — 
as the driving frequency was taken to be the same as the internal frequency of the in- 
finitely long continuous Josephson junctions. Yet, one needs to note that to simulate 
the governing equation numerically, the equation is discretized and solved on a finite 
interval, which implies that the system's internal frequency is likely different from the 
original equation. Therefore, p may not be necessarily zero. 

Treating p as a fitting parameter, we are able to find a good agreement between 
the numerics and the approximations for p « 0. Shown in the top panel of Fig. 16.21 
in dashed line is the approximation p.9p using p — 0.00607, where one can see that 
our approximation is in a good agreement as it is rather undistinguishable from the 
numerical result. In the middle panel of the same figure in dashed and dash-dotted 
line are the approximations for the driving amplitude h = 0.003 with p = 0.006 and 
p = 0.00665, respectively. The two values of p give a good approximation in different 
time intervals. It is surprising to see that the amplitude equation is still able to 
quantitatively capture the numerical result considering the large amplitude produced 
by the forcing, which is rather beyond the smallness assumption of the oscillation 
amplitude. 

In the bottom panel of Fig. 16.21 we plot the amplitude of the breathing mode in 
the — K junction case with h — 0.01, where one can see that the envelope of the 
oscillation amplitude tends to a constant. The dashed curve depicts our approximation 
from (|5.10p with p = —0.0015, where a good agreement is obtained. 

Considering the panels in Fig. 16. 2[ we observe that the mode in the two junction 
types does not oscillate with an unbounded or growing amplitude. After a while, there 
is a balance of energy input into the breathing mode due to the external drive and 
the radiative damping. The regular oscillation of the mode in the top panel indicates 
that the junction voltage vanishes, even when the driving frequency is the same as 
the system's eigenfrequency. This raises the question as to whether the breathing 
mode of a junction with a phase-shift can be excited further by increasing the driving 
amplitude to switch the junction to a nonzero voltage. To answer this question, 
we have solely used numerical simulations of (|1.2p as it is beyond our perturbation 
analysis. 

In the top left and right panels of Fig. 16. 3[ we present the average voltage (jl.ip 
with T = 100 as a function of the external driving amplitude h for the case of — tt — 
and — K junctions, respectively. One can clearly see that in both cases, there is a 
minimum amplitude above which the junction has a large nonzero voltage. For the 
first and second junction, the critical amplitude is respectively h « 0.34 and h w 0.1. 
The time dynamics of the transition from the superconducting state V ~ to a 
resistive state y ^ is shown in the bottom panels of the same figure. 

Form the panels, it is important to note that apparently the switch from a su- 
perconducting to a resistive state is not caused by the breathing mode, but rather 
because of the continuous wave background emitted by the breathing mode. It shows 
that the continuous wave becomes modulationally unstable and as typical dynamics 
breathers are created which later interact and destroy the breathing mode. Hence, we 
conclude that a breathing mode in these cases cannot be excited to make the junction 
resistive by applying an external drive, even with a relatively large driving amplitude. 

7. Conclusions. We have considered a spatially inhomogeneous sine-Gordon 
equation with a time-periodic drive, modeling a microwave driven long Josephson 
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junction with phase-shifts. Using muftiple scale expansions, we have shown that in 
an infinitely long Josephson junction, an external drive cannot excite the defect mode 
of a junction, i.e. a breathing mode, to switch the junction into a resistive state. For 
a small drive amplitude, there will be an energy balance between the energy input 
given by the external drive and the energy output due to so-called radiative damping 
experienced by the mode. When the external drive amplitude is large enough, the 
junction can indeed switch to a resistive state. Yet, this can also be caused by a 
modulational instability of the continuous wave emitted by the oscillating mode. 

Despite the agreement with the experiments obtained herein, our analysis is based 
on a simplified model. It is then of interest to extend the study to the case of (dc) 
driven long, but finite Josephson junctions with phase-shifts, as experimentally used in 
[131114] . In microwave-driven finite junctions, the boundaries can be a major external 
drive (see, e.g., [lllllO]), which is not present in the study here. A constant (dc) bias 
current, which is mentioned to play an important role in the measurements reported 
in [131 . is also not included in our current paper, even though the results presented 
herein should still hold for small enough constant drive. These are currently being 
studied and will be reported elsewhere. Another open problem that will be studied is 
the interaction of multiple defect modes |31j in a Josephson junction with phase-shifts. 
This is experimentally relevant as the so-called zigzag junctions have been successfully 
fabricated [18) . 
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Appendix. Explicit expressions. 
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The functions k^j and k^j in the expression of Vj (|2.17p . j = 1,2,3, are given by 



(^a)=2c.tan-i v/T+^(e^"+'^''^^(l + c.2) + 3e«Xov/r^) 



e"a; 2 w^" - 3 + 2 + 7 



(^^^ sin(2VTT^Xo)(2^+3)(.2 + i)2 ^ ^^^^^ ^ ^ 
4v 1 + a;2 2 
2 cos3 {VTTT^Xo) sin (\/r~Kj2Xo) tan"! [^J j^) + I + oj^) 



(^24)= 



2vr 



2 cos" (Vl + t^2Xo) (Vl tan-1 [^j^) + 1 + u;^) 



e"-2^«^i^(2w6_3_^_2^2^7^4) (1+ w2)'e3"-4^o/T^ 



2 



/ vr^ 



+e"Xo Vl-w^tan-i J^Y^i— _ 7^4 _ 2 ^2 _^ 3 _ 2 ^6 j ^ 
The functions £;j(Xo) and i?j(Xo) in 4^' (|431)-(|431) (j 0,2) are given by 



e(^»+"")Coi 2(2 + VT^7J^)e(^°+"")(2v^T^+3) 

(A.7) £;o(Xo) ~ ^ ^ g2(Xo+.o) VT^(1 + e2(Xo+xo))5 

(1 + Vl — w2)e2 Vl-t^^(Xo+2;o)+A'o+rE,) 

Vl-a;2(i + e2(^o+^o))5 

2 (2 - Vl - a;2)e(^« + "=''K2 ^ g g(Xo+a:o)(5+2 /T=^) 

Vl~r7j2(i + g2(Xo+Xo))5 
2(1 + 2 VI _ ^2y(Xo+xo)(2 

^ (2 + g2(Xo+2;o))5 ' 

e(-^°-=^o)Co2 2 e^-Xo+xo){2 x/T^-7)(2 + Vl - w^) 



(A.8)£;o(Xo) 



2 4.g2(Xo-xo) Vl - ^2(1 + e2(-^o-xo))5 



{A.9)E2{Xo) 
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1 _|_ g2(Xo+a:(3) 



((2 Vl -a;2 + i)e3(^o+rr„) _^ (2 Vl - - l)e^(-^''+^°)) VT^{Xo+xo) 

(1 + g2(J>fo+:ro))5 



2(1 + e2(^o+i;o))5 
(1 - \/l - a;2)e(^o+^»)(9+2vT^^) 

2((l + e2(-^«+^o))5 ' 



(2c^2g2(Xo-.„) ^ 2c^2 + ^i_4^2_i) g-yra^(Xo-.„)^ 

(A-10^2(Xo) = A l^g2(Xo-.o) ^ 



^g9{Xo-xo) _ 2 e3 (-^0-2:0) _ g(Xo-a:o)~) g-2 \/l-w^(Xo-rE„) 

2(l + e2(^o-^a))5 
2g-(Xo-i;o)(2\/T^;r2"-7) ^ ^^-(Xo+xo)i2VT^-3)^l _ ^^2 

2((l + e2(^o-:ro))5 



2(1 + e2(^o-a;o))5 

Vl - a;2(e-(^o-2:o)(2yT^i:^-9) _^ 4 g-(Xo-2:o)(2 TT^I^-T)^ 
2(1 + e2(^o-2:o))5 ■ 

The functions (7i(Xo), gi{Xo), and A^, j = 1, . . . , 4, in 



^ ^ {{XQ + Xo)c.2 + 1)7T37J2 _ 1^2 ^ ^)g(.o+Xo)(yT3^+2) 

(-2 + {xQ + Xq)c^4 _ + Xq - |)^2)yY^^gV73^(xo+Jfo) 

2a;Vl - u;2(^ + l)(w - l)(i + e2(^o+Xo) + i e4(:ro+Xo)) 



(2 + (-1 + + Xo)l^2)(^2 „ ^))g7CT(xo+^o) 

"^2wVl-u;2(a; + l)(a; - l)(i + e2(^o+-^") + i e4(^o+Xo) 
w2((xo + Xo)lu^ -xo-Xo- i)yr^7j2e(-o+x„)(4+vT^) 

^ 2w\/]~7j2(a; + l)(a; - l)(i + e2(^o+Xo) + l e4(a;o+Xo)) 



((-XQ -Xq- 1)c^2 ^ ^ ^ 1 )e(.o+Xo)(4+7T3ZJ2) 
"^2a;vT^^(w + l)(u; - l)(i + e2(»^o+Xo) + i g4(:.o+Xo)) ' 

(A.12) giiXo) 



((1 - (.To - Xo)oj^)VT^^ +1^^ - l)e(^^^^-2)(:.o-^o) 



wVl-w2(i + e-2(^o-Xo) + 1 e-4(:ro-Xo)) 
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2cjyr^^(w + l)(cj - l)(i + e-2(^-o--fo) + 1 e-4(:.o-x„)) 

2wVl - W2(CJ + - l)(i + e-2(^0-X„) + 1 g-4(:ro-Xo)) 



(c^ + l)(-a;o + + + ^)e-(-^o+^o)(^^-4) 
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